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CHAPTER  I 
INTRODUCTION 


In  recent  years  many  new  techniques  have  been  developed  for 
analyzing  systems  which  could  not  be  adequately  handled  by  classical 
methods.  Two  of  the  most  powerful  of  these  theories  are  the  moment 
method  (MM)  and  the  geometrical  theory  of  diffraction  (GTD).  The 
moment  method  is  a  variational  technique  for  solving  integral  equa¬ 
tions  by  reducing  them  to  a  system  of  N  algebraic  equations  with 
N  unknown  coefficients.  Since  most  electromagnetic  scattering  prob¬ 
lems  can  be  formulated  as  an  integral  equation,  the  moment  method 
finds  wide  application.  Due  to  the  fact  that  the  number  of  unknowns 
is  more  or  less  directly  related  to  the  size  of  the  scattering  sys¬ 
tem,  the  moment  method  is  usually  limited  to  systems  which  are  not 
electrically  large. 

The  GTD  is  an  asymptotic  technique  which  takes  advantage  of 
the  fact  that  electromagnetic  radiation  behaves  as  a  ray  optical 
phenomenon.  The  usefulness  of  GTD  rests  on  the  idea  that  interaction 
between  rays  and  objects  are  local  occur ances  and  thus  it  is  not 
necessary  to  integrate  currents  over  an  entire  structure  as  is  done 
In  classical  techniques.  This  local  nature  of  the  interaction, 
allows  complex  geometries  to  be  built  from  canonical  forms  such  as 
flat  plates,  cylinders,  wedges,  cones,  etc.  However,  because  of 
its  asymptotic  formulation,  GTD  is  primarily  useful  for  systems  which 
are  electrically  large.  Thus  there  exists  a  gap  in  the  tools  avail¬ 
able  for  analyzing  scattering  systems.  Systems  which  contain  an 
electrically  small  scatterer  in  the  presence  of  a  large  body  cannot 
be  adequately  analyzed  by  either  the  moment  method  or  GTD 

The  hybrid  method  presented  in  this  paper  presents  a  technique 
for  analyzing  such  a  system.  The  work,  presented  herein  is  based 
on  previous  papers  by  Thiele  and  Newhouse  [7]  and  by  Ekelman  Cl], 

The  approach  taken  is  to  cast  the  small  scatterer  (built  of  thin 
wires)  in  the  moment  method  format  of  Richmond  L8]  and  then  add  per¬ 
turbation  terms  due  to  interaction  with  the  nearby  large  body  using 
GTD.  Thus  the  solution  is  basically  a  GTD  perturbation  of  the  moment 
method.  The  particular  form  of  the  GTD  used  in  this  work  is  a  new 
ansatz  called  the  uniform  theory  of  diffraction  (UTD)  reported  by 
Pathak  et  al  [3].  The  UTD  avoids  previous  difficulties  encountered 
with  GTD  when  a  ray  hits  a  curved  edge  at  grazing  incidence. 
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The  hybrid  method  is  a  general  technique  which  is  applicable 
to  a  large  class  of  problems.  The  particular  solutions  presented 
herein  extend  the  initial  work  done  by  Ekelman,  providing  further 
examples  of  the  utility  of  the  hybrid  technique  and  greatly  broaden 
the  scope  of  problems  for  which  solutions  are  currently  available. 
The  first  solution  given  is  for  two  piecewise  sinusoidal  dipoles 
arbitrarily  positioned  around  an  infinite  circular  cylinder.  This, 
along  with  Richmond's  work,  provides  a  technique  for  analyzing  any 
arbitrary  combination  of  thin  wire  antennas  near  a  conducting  cyl¬ 
inder.  The  second  solution  provides  a  method  for  analyzing  antennas 
mounted  on  the  cylinder,  such  as  a  monopole.  This  solution  is  ob¬ 
tained  as  a  perturbation  of  image  theory.  This  monopole  solution 
is,  in  principle,  as  general  as  the  rest  of  the  hybrid  theory,  how¬ 
ever  the  GTD  expression  needed  for  a  completely  general  monopole 
analysis  were  not  available.  Thus  the  data  presented  does  not  in¬ 
clude  the  full  GTD. 

Chapter  II  contains  the  development  of  the  moment  method  as 
it  is  used  in  this  work.  Chapter  III  presents  the  essential  points 
of  the  GTD  (UTD)  and  included  the  UTD  expressions  of  Pathak  et  al. 
The  hybrid  method  is  covered  in  detail  in  Chapter  IV.  Results  of 
the  two  dipole  problem  are  presented  in  Chapter  V  and  Chapter  VI 
contains  the  theory  and  implementation  notes  on  the  monopole  problem 
from  conclusing  remarks  are  contained  in  Cahpter  VII. 


CHAPTER  II 
THE  MOMENT  METHOD 


The  moment  method  is  a  procedure  for  solving  an  integral  equa¬ 
tion  of  the  first  kind  as  in  Equation  (1) 

f I(r')K(r,r')dr'  =  -E1 (r)  (1) 

over  a  structure 


Almost  all  electromagnetic  scattering  problems  can  be  put  into 
the  form  of  Equation  (1)  and  then,  in  principle,  solved  by  moment- 
method  techniques.  Since  the  moment  method  is  such  an  extremely 
powerful  mathematical  tool,  a  completely  general  discussion  would 
be  quite  abstract  and  thus  not  entirely  suitable  for  the  purposes 
of  this  paper.  We  will  therefore  consider  a  specialized  case  of 
the  moment  method.  We  will  consider  an  arbitrary  structure  built 
from  electrically  thin  wires  (wire  radius  <  .006a).  Since  a  large 
portion  of  antennas  in  common  use  today  are  made  of  or  can  be  modeled 
by  thin  wires,  this  solution  will  be  applicable  to  a  large  class 
of  problems.  For  such  a  structure,  the  scattering  is  given  by 
Pocklington's  integral  equation  (2)  for  a  z  directed  segment 

/I(z')[v2G(z,z')  +  02G(z,z')]dz'  =  -Ej(z)  (2) 

0  over  structure 

where  G(z.z')  is  the  free  space  Green's  function 
e-je| z-z'| 

G(z,z')  ■  -  )z-z.|  {3) 

The  function  I ( z * )  is  the  current  in  the  wires  and  Ejz)  is 
the  incident  or  impressed  field  at  the  surface  of  the  wire. 

Note  that  the  volume  integral  of  (1)  has  become  a  line  integral 
in  (2)  since  the  current  on  a  thin  wire  is  assumed  to  be  a  line 
source.  A  typical  example  of  the  type  of  problem  which  can  be  solved 
by  (2)  is  shown  in  Figure  1  where  a  thin  wire  h  wavelengths  long 
(h«l)  is  illuminated  by  a  plane  wave  Thp  quantity  E’(z)  is  the 
electric  field  incident  on  the  wire  and  I(z)  is  the  anxnown  currents 
induced  by  the  incident  field. 
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The  simple  example  in  Figure  1  can  be  solved  by  assuming  a  cur¬ 
rent 


I(z) 


I0sinn(^  -|z|  ) 
sin(^  ) 


(4) 


Figure  1.  Thin  wire  illuminated  by  a  plane  wave. 


Substituting  (4)  into  (2)  along  with  the  known  function  E’(z) 
will  give  I  and  thus  the  problem  is  solved.  For  a  more  complex 
structure,  the  form  of  the  current  is  not  known  and  thus  no  assump 
tion  can  be  made  as  in  (4),  in  this  situation  the  moment  method  is 
used. 


The  moment-method  technique  produces  numerical  solution  to  an 
equation  such  as  (2)  by  transforming  the  integral  equation  into  a 
set  of  linear  algebraic  equations  as  follows.  Consider  an  arbitrary, 
complex  thin  wire  structure.  The  unknown  currents  on  such  a  struc¬ 
ture  can  be  expanded  as  a  series  of  expansion  functions 

N 

I(z)  =  l  InJ  (z)  (5) 

n=l  "  n 


where  the  expansion  functions  J  ( z)  are  arbitrarily  chosen  (as  in 
the  example  for  Figure  1)  and  defined  over  a  particular  portion 
of  the  structure,  and  zero  elsewhere.  These  functions  along  with 
the  weighting  coefficients  I  ,  describe  an  arbitrary  current 
distribution  as  N-*».  ObviouSly,  for  numerical  computation  it  is 
desireable  to  make  N  as  small  as  possible  and  still  obtain  accurate 
results.  The  choice  of  J  will  be  discussed  shortly. 
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Substituting  (5)  into  (1)  gives 


N  i 

InJn(z')KU,z')dz'  =-E’(z)  (6) 

n=l 

Now  all  of  the  terms  inside  the  integral  are  known,  and  thus  the 
integral  can  be  evaluated  to  give 

X  W2>  *  -<<z>  <7> 

where 

Fn(z)  Jn(z)K(2,2’)dzt  (8) 


Now  the  integral  has  been  removed  but  the  result  is  one  equation 
with  (N)  unknown.  The  (N)  equation  needed  to  solve  (7)  and  be  pro¬ 
duced  using  weighted  residuals. 

Weighted  Residuals 

The  well  known  boundary  conditions  at  the  surface  of  a  perfectly 
conducting  wire  are 


where  Ej  and  e!:  are  respectively  the  tangential  incident  and  scat¬ 
tered  fields.  Since  the  assumed  form  of  the  current  (Equation  (5)) 
is  in  general  not  exact  unless  N=<*>,  then  the  residual  defined  in 
Equation  (10)  will  not  be  zero 

R(z)  =  Ej(z)  +  Ej(z)  (10) 

for  all  z. 

In  order  to  preserve  the  boundary  condition  as  nearly  as  pos¬ 
sible,  the  residual  is  forced  to  zero  in  an  average  sense.  It  is 
hoped  that  In  so  doing,  the  result  will  be  as  near  to  the  exact  an¬ 
swer  as  possible.  Using  a  weighted  average  gives. 
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JW  (z)R(z)dz  =  0  J  n=l,2,3. . 
over  structure 


(ID 


here  W  (z)  are  arbitrary  weighting  functions.  The  choice  of  these 
will  b?  discussed  shortly. 

Substituting  (11)  and  (10)  into  (7)  gives 


over 


■fVz)  l  InFn(z)dz  =  f  E](z)W  (z)dz 
structure  n=l  over  structure 


(12) 


m  =  1,2, ...N 

which  will  yield  the  required  N  equations.  Equation  (12)  can  be 
written  in  a  simpler  form  by  substituting 


Vm  =  '  Lr  stru!ture*<Z>dZ 


Zmn  a  1  Vz)Fn(2)<lz 
over  structure 


Then 


N 

V  I  Z  =  V  m=l,2,...N 
n  mn  m  ’  ’ 


(13) 

(14) 


(15) 


or  in  matrix  notation 
[V]  =  [Z][I  ] 


(16) 


Obviously  (16)  can  be  solved  by  suitable  techniques,  thus  yield¬ 
ing  the  required  coefficients  ( In) . 

The  choice  of  J  ( z)  and  W (z)  have  a  great  deal  to  do  with  the 
number  of  terms  (N)  needed  to  get  an  accurate  result.  In  general  the 
more  closely  these  functions  model  the  actual  current  distribution 
on  the  wire,  the  smaller  (N)  can  be.  When  J_(z)  and  W  (z)  are  the 
same,  one  obtains  a  special  case  of  the  moment  method  called  Galerkin's 
method. 
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The  Reaction  Integral 

The  preceeding  development  of  the  moment  method,  while  mathe¬ 
matically  correct,  gives  little  physical  insight  into  the  quantities 
Z  and  V  (Equations (13),  (14)).  These  quantities  can  be  shown 
tonhave  physical  significance  as  follows. 

In  1954,  Rumseytll]  introduced  a  physical  observable  called 
reaction,  as 

ff  ( J*E  -M*H  )ds  -  rr  J-E.-M  -H.)ds  =  0  (17) 

J  J  m  s  m  s  J^mimv  '  ' 

where  (E  ,H  )  denotes  the  field  generated  by  the  equivalent  surface 
current  densities  (J  M  )  of  a  test  source  (J  M  ),  as  in  figure 
2.  Also  (E.,H.)  are  the  rields  of  an  externa  i  source  (J1  M1).  Phys¬ 
ically,  Equation  (17)  states  that  there  is  zero  reaction  or  coupling 
between  the  test  source  and  the  sum  of  the  incident  and  scattered 
fields.  This  is  equivalent  to  the  weighted  residual  interpretation 
of  Equation  (11).  For  the  case  of  thin  wires  and  electric  currents, 

M =  0  and  the  surface  integrals  reduce  to  line  integrals,  thus 
(17)  becomes 

/(Vrs)dz  =  /  ‘  ri)ds  (18) 

which  is  the  same  form  as  (12).  Thus  the  weight  function  can  be 
interpreted  as  test  currents  or  test  sources.  If  the  expansion  func¬ 
tions  are  taken  to  be  at  the  center  of  the  perfectly  conducting 
wires,  and  the  surface  s  in  Figure  2  is  taken  to  be  the  surface  of 
the  wire,  then  Z  is  merely  the  coupling  between  the  mth  test  func¬ 
tion  and  the  nth  expansion  function. 


Figure  2.  Test  source  (IT  m  )  inside  conducting 
.  ,  ,  _m  m  J 

body  surface  S. 
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Piecewise  Sinusoidal  Galerkin  Method 


The  piecewise  sinusoid.  Figure  3,  used  by  Richmond  [8  ]is  one 
of  the  most  useful  basis  functions  for  thin  wires.  It  may  be  ex¬ 
pressed  mathematically  for  each  z-directed  segment  by 


I  ( z ) 


I^sinhy^-zJ  +  ^sinhyU-z^) 
sTnyd- 


(19) 


where  I,  and  I„  are  the  endpoint  currents,  y  is  the  complex  propa¬ 
gation  constant  of  the  medium,  and  d=z2-z,  is  the  source  length. 
Piecewise  sinusoids  are  efficient  since  very  few  are  needed  to  ap¬ 
proximate  the  expected  current  distribution  as  in  Figure  4. 


k 


Figure  3.  Typical  two  segment  dipole  with 

piecewise  sinusoidal  distribution. 


Figure  4.  Current  distribution  modeled  by  overlapping 
piecewise  sinusoidal  modes. 
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Note  that  each  piecewise  sinusoid  function  spans  two  connected 
segments  and  that  each  interior  segment  contains  two  piecewise  sinu¬ 
soidal  functions.  Thus  it  is  convenient  to  consider  each  piecewise 
sinusoid  as  a  dipole  with  its  source  at  the  junction  between  the 
segments.  Consider  the  situation  wherein  an  arbitrary  number  of 
segments  are  coincident  with  the  z-axis.  For  example,  for  an  antenna 
with  N  segments,  the  mn-th  impedance  matrix  element  is 


?  <2-Vi> 

i,  s1nIVrzJ 


Z-E'dZ 


zm+l  sin  (z.,,-z) 

c 

The  field  T  needed  to  evaluate  Equation  (20)  may  be  found  in  a 
straightforward  manner,  and  the  derivation  is  included  in  many  electro¬ 
magnetic  texts  including L  9].  The  geometry  is  shown  in  Figure  5. 

For  the  piecewise  sinusoidal  excitation  current  Es  is 


4-rrpS 


\(  '^1 

Vnh^dl^1!6 


"11*1  -Yr2> 

-I2e  2J 


inhyd 


+  (IjCosh  I-I2)e  1cos61 

-Yr?  I 
+  ( I^cosh  d-I^e  cosfU 


Ez  =  4  slnhyd  [(IrI2cosh  d>  TP 


+  (I2-IjC0sh  d) 

where  n  is  the  impedance  of  the  medium.  The  evaluation  of  Equation 
(20)  may  be  carried  out  without  difficulty  by  numerical  integration 
or  by  Si  and  Ci  integrals  for  some  special  geometries. 


The  piecewise  sinusoidal  Galerkin  method  presented  here  is  the 
procedure  used  for  the  th in-wire  antennas  of  this  paper.  The  im¬ 
portance  of  this  particular  specialization  of  the  moment  method  for 
the  hybrid  technique  will  be  discussed  in  Chapter  IV. 
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CHAPTER  III 

GEOMETRICAL  THEORY  OF  DIFFRACTION 


While  the  moment  method  is  a  technique  primarily  for  low  fre¬ 
quency,  the  geometrical  theory  of  diffraction  (GTD)  is  useful  for 
high  frequency,  or  situations  where  the  scatterer  is  electrirally 
large.  The  foundation  of  GTD  rests  on  two  basic  principles. 

i.  Electromagnetic  propagation  is  a  ray  type  phenomenon,  that 
is,  a  wave  travelling  from  point  a  in  space  to  point  b, 
will  follow  a  geodesic  path.  In  general,  there  is  more 
than  one  such  path  and  the  contribution  from  all  such  paths 
must  be  considered  before  an  accurate  solution  can  be 
reached.  Several  examples  are  shown  in  Figure  6. 

ii.  Secondly,  all  interaction  between  rays  and  objects  are 
localized  phenomenon.  Thus  it  is  not  necessary  to  inte¬ 
grate  the  currents  over  an  entire  structure  in  order  to 
obtain  the  scattered  field.  It  is  only  necessary  to  de¬ 
termine  the  incident  ray  and  the  geometry  of  the  object 
at  the  points  where  the  ray  is  in  "contact"  with  the  sur¬ 
face. 

The  localized  nature  of  the  interaction  and  the  ray  optic  scat¬ 
tering  at  once  enable  one  to  model  a  complex  structure  from  simple 
geometric  shapes,  such  as  cylinders  and  plates.  Once  this  model 
has  been  constructed  one  can  quickly  evaluate  the  components  of  the 
structure  which  will  have  a  significant  effect  on  the  scattered 
field.  For  most  structures  there  are  an  infinite  number  of  geodesic 
ray  paths  from  point  a  to  point  b,  however  the  nature  of  the  dif¬ 
fraction  allows  one  to  quickly  determine  the  significant  rays  in 
order  to  arrive  at  an  engineering  solution.  Assuming  the  scatterer 
does  not  hide  the  source,  the  dominant  field  is  produced  by  the  di¬ 
rect  ray  from  the  source  to  the  observer.  The  next  most  important 
term  is  the  ray  reflected  from  the  scatterer.  This  field  is  termed 
the  geometrical  optics  field. 

A.  The  Geometrical  Optics  Field 

Although  the  GTD  can  handle  a  general  astigmatic  wave  front, 
only  spherical  waves  will  be  discussed  here  for  reasons  that  will 
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be  apparent  later.  The  geometry,  shown  in  Figure  7,  consists  of 
an  arbitrarily  curved  smooth,  perfectly  conducting  surface  with  near¬ 
by  source  and  observation  points.  The  ray  path  is  from  the  source 
point  to  the  reflection  point,  and  then  to  the  observation  point. 

The  reflected  field  Er(s)  of  a  curved  surface  S  is 


Er(s)  -  E\QR) 


r  r 

plp2 _ 

(Pj+s)  (p^5) 


(23) 


where  E^Q-)  is  the  incident  field  at  the  reflection  point  QR.  In 
(23),  as  in  all  of  the  formulas  given,  the  eJu)t  time  dependence 
is  assumed  but  not  shown. 


OBSERVATION 


R  is  the  dyadic  reflection  coefficient  expressed  as 


%  *1  Af*  Al  At' 

H  *  e  e  -  e  e 


(24) 


The  unit  vectors  e^and.e^are  parallel  to  the  plane  of  incidence 
and  the  unit  vectors  e?,  ef  are  perpendicular  to  the  plane  of  inci¬ 
dence.  The  reflection  point  Q  must  be  determined  from  the  law  of 
reflection  r 


-f 


n  =  s 


(25) 


and 


txn=sxn  (26) 

or  simply  stated,  (25)  states  that  the  angle  of  incidence  is  equal 
to  the  angle  of  reflection,  and  (26)  states  that  the  incident,  re¬ 
flected  and  surface  normal  vector  all  lie  in  the  same  plane.  The 
quantities  p.p„  are  the  principal  radii  of  curvature  of  the  reflected 
wave.  These^ttere  given  by  Kouyoumjian  and  Pathak  [4]  for  spherical 
wavefronts  as 


1 

pl,2 


s 1 


cose 


sin20. 


2  2  2 
sin  e2  since1'j  4 

+  ~Tq~  +  -T— ] 


(27) 


where  s'  is  the  distance  from  the  source  to  Q  ,  0  is  the  angle  be¬ 
tween  the  direction  of  the  incident  ray  I  andrunit  vector  e,  and 
6 2  is  the  angle  between  incident  ray  I  and  unit  vecotr  e^. 

B.  Curved  Surface  Diffraction  and  UTD 


The  geometrical  optics  field  is  valid  only  when  the  source  and 
observation  points  are  not  hidden  by  the  structure  and  the  observa¬ 
tion  point  is  well  away  from  the  shadow  boundary.  A  more  general 
case  is  when  the  observation  point  can  be  anywhere  in  the  space  around 
the  structure.  This  problem  is  much  more  complex.  Although  the 
results  given  in  this  section  are  valid  for  any  convex  surface  we 
will  limit  our  discussion  to  an  elliptical  cylinder  which  is  infinite 
in  extent.  This  will  serve  to  demonstrate  the  essential  points  of 
the  solution  without  adding  unnecessary  complication.  The  regions 
of  interest  around  an  elliptic  cylinder  are  shown  in  Figure  8. 

The  uniform  GTD  (UTD)  ansatz  of  Pathak  et  al.  [3]  was  chosen  because 
of  its  good  continuity  through  the  shadow  boundary  where  previous 
theories  failed.  This  solution  is  valid  when  the  source  and/or  ob¬ 
servation  points  are  in  regions  I,  II,  and  III. 
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PL  ■  FIELD  POINT  IN  LIT  ZONE 
Ps>  FIELD  POINT  IN  SHADOW  ZONE 


Figure  8.  The  rays  and  regions  associated  with  scattering 
by  a  smooth  convex  cylinder. 


The  geometries  to  be  considered  are  shown  in  Figures  7  and  9. 
Figure  7  shows  the  geometry  to  be  considered  when  the  observation 
point  is  in  the  lit  zone  (reflection).  Figure  9  shows  the  geometry 
to  be  considered  when  the  observation  point  is  in  the  shadow  region 
(creeping  wave  diffraction).  Note  that  for  this  case  the  ray  couples 
to  the  surface  at  point  Q.  and  then  "creeps"  along  the  curved  sur¬ 
face  to  the  shed  point  Q_.  For  the  reflection  problem  the  scattered 
field  is  computed  as  in  the  geometrical  optics  case  but  with  a  gen¬ 
eralized  E  dyadic.  In  the  creeping  wave  problem,  the  scattered  field 
is  computed  from  the  incident  field  and  the  F  dyadic.  Referring 
to  the  appropriate  geometry,  the  1?  and  T  dyadics  are  defined  as 

=  Rc@1er  +  R.  eV  (28) 

s  i  i  n  ii  H 


T  ■  TA62  ♦  Wj 


(29) 


where 
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Figure  9.  Diffraction  by  a  smooth  curved  surface. 
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shadow  region. 


***  . 


with 


=-2m(QR)  f(QR)  ”1/3cos61 ; 


,(Qo) 


f(Q«)=1+  £?X)  cos2ei;  !  dt'  ^  ; 

t  K 


m 


m(Q) 


kp  (Q) 


1/3 


*2 

t  -  /  df  ; 
*1 


(32a, b,c) 
(32d,e) 


XL  =  2k lL  cos2©1  ; 


xd  -  k  lV)2 

2m(Q1)m(Q2) 


(32f,g) 


Thus  with  the  "R  and  T  dyadics  computed,  the  scattered  fields  can 
be  computed  by 


e(pl) 


r 


(PL)+r1(QR)*I 


r  r 

plp2 _ 

(pK)(p2+sr) 


for  P,  in 
the  jrt  region 


E(PS)  .  Ei(Q1).T 


d 

p2 


S-Jks 


(P2+s°) 


(33) 


for  P  in  the 
shadoft  region 


(34) 


Note  that  in  the  above  equation 


[l-F(a)]  -*•  0  in  the  region  exterior  to  zone  II 


PAS)  ~  MR 


lps  . 

re  >n 


zone  I 


(35) 


Thus  the  above  solution  reduces  to  the  geometrical  optics  solution 
in  the  visible  region  away  from  the  shadow  boundary  (zone  I). 
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The  parameters  LL  and  Ld  are  distance  parameters  relating  to 
the  relative  sizes  of  s,s'.  for  a  spherical  wave  they  are  given 
by 


L  _  ss 

d  _  ss 
s+s 


s+s 1 


While  Equations  (33)  and  (34)  are  valid  and  continuous  at  the  ^hadow 
boundary  they  are  only  true  in  the  limit  since  the  parameters  r, 
and  C  are  both  zero  on  the  shadow  boundary.  A  limit  form  of  the 
solution  which  can  be  employed  in  this  region  is 


‘  Ei(PSB) 


E1 (Qj ) 


JTL 

jr  +  >"(0^ 


•Ei(Q1) 


J7 

T 


+  m(Q1) 


The  Ps(£)  functions  of  Equations  (30,31) 


the  related  functions  p*(?)  and  q*(£)  by 

fp*(6) 


Ps(«)  =  i 

h 


q*(6) 


r  e 


-i  H 
J  4 


-j  JL 
.  3  4 


26  JW 


are  defined  in  terms  of 


(37) 


where  the  p( 6)  and  q(6)  functions,  defined  by  Logan  [6]  are  plotted 
in  Figures  10a, b  over  the  region  -3<o<  2.  More  complete  descriptions 
of  these  functions  are  given  in  [2]  and  [5].  The  transition  function 
F ( <s )  which  also  appears  in  formulas  (30),  (31)  is  defined  [4]  in 
terms  of  the  Fresnel  integral  as 


F( 6)  =  2j  6  ejx 


dx  . 


(38) 


i 
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CHAPTER  IV 
THE  HYBRID  METHOD 


The  hybrid  method  to  be  presented  here  is  a  technique  for  an¬ 
alyzing  thin  wire  antennas  which  are  not  electrically  large  in  the 
presence  of  a  body  or  scatterer  which  is  too  large  to  be  handled 
by  the  moment  method.  This  is  accomplished  by  computing  a  second 
moment  method  impedance  matrix  proposed  by  Thiele  et  al  [7],  This 
second  matrix  Ia Z J  is  a  perturbation  to  the  original  [Z]  and  includes 
the  coupling  between  the  test  and  expansion  modes  due  to  rays  re¬ 
flected  and  diffracted  from  the  nearby  structure. 


Using  the  GTD  ray  optics  technique,  the  electric  field  incident 
on  a  wire  scatterer  due  to  radiation  from  another  wire  scatterer, 
can  be  written  as 


(39) 


where  E°  in. the  incident  field  due  to  the  direct  path  (free  space 
case)  and  E denote  the  fields  due  to  the  various  rays  which 
are  reflected  orsdiffracted  by  the  nearby  body.  Note  that  fog  the 
case  where  the  body  hides  the  source  from  the  observer  then  E°  = 

0.  Using  Equation  (39)  the  component  Z  of  the  [Z  ]matrix  cin  be 
written  mn 

Zmn  =  tfy)  +  Es(£)  +  Es  **’  Es  ]  *  (40) 


which  can  be  written  as 

Z  =  Z '  +  &Z 

mn  mn  mn 


(41) 


where 


Z 


mn 


=  /  E°(d)  •  I(*)dH 


(42) 


and 

^mn  =  +  Es<*>  *•*  E^>  ’  IW*  («) 


and  thus  in  matrix  form 


[Z]  =  [Z1]  +[AZ]  (44) 

Where  [ Z ' ]  is  computed  using  the  direct  (free  space)  path  between 
source  and  observation  point  and  [AZ  ]  is  computed  using  the  ray  paths 
which  interact  with  the  large  body.  Thus  the  effect  of  the  large 
body  are  included  in  the  moment  formulation  and  the  solution  obtained 
will  therefore  contain  these  interactions.  One  problem  which  arises 
from  Equation  (39)  is  that  the  GTD  requires  that  the  field  generated 
by  the  sources  contain  no  ray  path  components.  In  general  this  con¬ 
dition  cannot  be  met.  Fortunately,  the  PWS  expansion  function  used 
by  Richmond  can  be  shown  to  have  only  transverse  field  components 
as  follows; 


Eckelman  [  1  ]  observed  that  the  fields  of  a  piecewise  sinusoidal 
test  dipole  can  be  broken  up  into  contributions  from  the  endpoints 
and  the  driving  points.  These  fields  have  spherical  wavefronts  and 
thus  are  in  the  required  form  for  using  them  as  sources  with  the 
GTD  (UTD).  With  reference  to  Figure  5,  the  exact  field  of  a  PWS 
dipole  is 


E 


1 


E 


2 


4  7rsin(kd) 


-jkR-i 

e 


jnI02cos(kd)  e‘jkR2 
2  -n  s  f  i nCkdT  sT n  a,~  R~  ~ 


(45a) 


(45b) 


f  -~Jn  lo 
j  4ti  sfinTkd) 


1 

sinQj 


(45c) 


Thus  using  the  above  result,  the  coupling  between  a  test  source  and 
an  expansion  mode  proceeds  as  follows: 


A  term  Z  of  the  Z  matrix  is  computed  in  the  standard  form 
using  formulamv20) ,  this  gives  the  coupling  due  to  the  direct  ray 
path.  Note  that  if  the  test  and  expansion  modes  are  hidden  by  the 
cylinder  then  this  Z  will  be  zero.  The  elements  of  £  are  com¬ 
puted  also  using  formula  (20)  where  E^  is  the  field  scattered  from 
the  cylinder.  This  includes  all  significant  fields  from  all  ray 
paths  of  interest.  For  this  paper,  these  ray  paths  are  defined  by 
the  statements: 

i)  if  the  observation  point  is  not  shadowed,  only  the  reflected 
field  is  included. 
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ii)  if  the  observation  point  is  shadowed  the  rays  corresponding 
to  creeping  around  the  cylinders  in  a  clockwise  and  a 
counter-clockwise  sense. 

These  are  the  significant  terms  for  an  infinite  cylinder.  If  the 
cylinder  were  not  infinite,  scattering  from  the  end  caps  would  have 
to  be  included. 

The  total  impedance  matrix  is  then  computed  as 
C  Z]  -  [Z‘  ]  +  [a Z] 

Once  the  [V]  matrix  has  been  computed  (in  the  usual  manner),  then 
[I]  can  be  found  by  matrix  inversion  or  any  other  suitable  technique. 
Once  again,  it  is  important  to  note  the  significance  of  using 
Richmond's  piecewise  sinusoidal  test  function  for  the  test  sources. 

As  a  result  of  the  endpoint  separation  technique  of  Eckelman  each 
term  in[aZ]can  be  computed  by  integrating  only  over  the  expansion 
mode,  and  using  the  field  radiated  by  the  endpoints  of  the  test  func¬ 
tion.  For  a  general  test  function,  two  problems  would  arise  which 
would  greatly  reduce  the  numerical  efficiency  and  accuracy  of  the 
process. 

1)  The  expression  for  the  coupling  would  contain  a  double 
integral:  one  over  the  source,  to  compute  the  incident 
field  at  a  point  on  the  observer,  and  one  over  the  observer 
to  compute  the  actual  coupling. 

2)  The  fields  radiated  by  the  source  would  not  necessarily 

be  transverse  to  the  direction  of  propagation,  w'-ich  cannot 
be  handled  using  GTD. 


CHAPTER  V 

DIPOLES  NEAR  AN  ELLIPTIC  CYLINDER 


The  first  application  of  the  theory  presented  in  Chapters  II, 
III,  and  IV  is  the  analysis  of  dipoles  near  an  infinite,  perfectly 
conducting,  elliptic  cylinder.  The  elliptic  cylinder  was  chosen 
because  it  has  found  wide  application  in  the  past.  It  can  be  used 
to  model  aircraft  fuselages  and  engines  as  well  as  antenna  masts 
on  ships  and  land  based  systems. 

The  dipole  is  of  prime  interest  because  as  shown  in  the  earlier 
chapter  on  the  moment  method,  the  [Z]  matrix  can  be  obtained  by  com¬ 
puting  the  coupling  between  the  mth  test  function  and  the  nth  ex¬ 
pansion  mode.  In  the  case  of  Richmond's  thin  wire  code,  these  modes 
are  dipoles  with  a  piecewise  sinusoidal  current  distribution.  Thus 
once  one  can  calculate  the  near  field  coupling  between  dipoles  near 
an  object,  any  thin  wire  antenna  near  the  object  can  be  analyzed. 

A.  Overview 


The  geometry  to  be  considered  consists  of  two  piecewise  sinu¬ 
soidal  dipoles  located  in  the  vicinity  of  an  infinite,  perfectly 
conducting  elliptic  cylinder  which  is  aligned  with  the  z  axis.  One 
dipole  will  be  referred  to  as  the  source;  the  other  dipole  will  be 
called  the  observer.  These  designations  are  arbitrary  since  by  reci¬ 
procity  the  source  and  observers  can  be  interchanged  without  affecting 
the  coupling.  Coupling  will  be  computed  using  Equation  (20).  The 
integration  over  the  observer  will  be  carried  out  numerically  by 
dividing  the  observer  up  into  subsection  and  using  Simpson's  rule 
integration.  The  fields  radiated  by  the  source  will  be  spherical 
waves  eminating  from  the  driving  point  and  the  two  endpoints.  For 
rays  which  interact  with  the  cylinder,  the  diffraction  or  reflection 
coefficients  will  be  computed  using  the  UTD  presented  in  Chapter 
III.  Examination  of  the  UTD  expressions  shows  that  it  is  first 
necessary  to  determine  whether  or  not  the  cylinder  hides  the  source 
from  the  observer. 

B.  Location  of  Observer 


The  question  of  determining  whether  the  cylinder  hides  the 
source  can  be  resolved  by  determining  whether  or  not  a  direct  ray 
from  source  to  observer  will  intersect  the  cylinder.  Since  the  cy¬ 
linder  is  infinite  in  length,  the  geometry  can  be  reduced  to  a  two 
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dimensional  problem  in  the  x,  y  plane.  Consider  the  geometry  shown 
in  Figure  11.  The  equation  of  the  line  (xsys)(xcyc)  can  be  written 


y  =  mx  +  c 


(46) 


where 


m  = 


V4 

VXr 


and 


c  -  ys-mx$ 


(47a, b) 


the  equation  for  an  ellipse  is 
2 

4  =  1  (481 

b^ 

thus  Equations  (46),  (47a, b)  and  (48)  can  be  combined  to  give 

(491 

then  using  the  quadratic  formul* 


x2  -2 
a2  b‘ 


m2  +  4  x2  +  2mcx  +  (c2-b2l  =  0 
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2  2  I?  +^JlcC-*C) 


2  .2, 


-2mc  -  / 4  me  - 


'1,2 


^m2  + 
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(50) 


If  x,  _  are  complex  conjugates,  then  no  solution  exists  in  the  x,y 
plane’and  thus  the  line  does  not  intersect  the  cylinder  at  all.  If 
x,  „  exist  then  the  situation  shown  in  Figure  12  can  be  detected 
as’rollows:  If  x,  <  x,  ?  <  x  or  x  <  x,  ?  <  x$  then  the  source 
is  hidden,  otherwise  aic6nditTon  1 i Ke  Figure  12  exists  and  the  source 
is  obviously  not  hidden. 

C.  Spectral  Point  Calculation 

If  the  wave  is  not  hidden,  then  it  is  necessary,  to  find  the 
reflection  point.  MarhefkaC  2]  outlined  a  procedure  for  finding 
the  reflection  point  on  an  elliptic  cylinder  assuming  the  source 
point  and  the  scatter  direction  were  known.  Ekelman  modified  this 
procedure  for  the  case  where  the  source  point  and  observation  point 
were  known.  This  procedure,  while  accurate,  is  inefficient.  The 
solution  presented  here  works  well  for  a  large  class  of  problems 
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r  ~ 


when  neither  the  source  nor  observation  point  are  in  the  far  field 
of  the  cylinder. 

The  geometry  of  the  problem  is  shown  in  Figure  13.  The  source 
and  observation  points  are  in  the  vicinity  of  an  elliptic  cylinder 
and  the  reflection  point  (x  y  z  )  on  the  surface  of  the  cylinder 
is  desired,  and  it  is  assumedrtnat  the  cylinder  does  not  hide  the 
source  from  the  observation  point.  Marhefka  [2]  formulates  the  prob¬ 
lem  using  the  laws  of  reflection  and  shows  that  the  problem  can  be 
decomposed  into  a  two  dimensional  problem  in  x  and  y,  and  a  one  dimen¬ 
sional  problem  in  z. 

Consider  the  two  dimensional  problems  shown  in  Figure  14. 

The  incident  and  observation  vectors  can  be  written 

T(v,z)  =  Ix(v)x  +  Iy(v)y  +  I2(z)z 

=  (acosv-xs)x  +  (bsinv-y$)y  +  (z-zs)z  (51) 

and 

d(v,z)  =  dx(v)x  +  dy(v)y  +  dz(z)z 

=  (x  -acosv)x  +  (y  -bsinv)y  +  (z  -z)z  .  (52) 

C  u  L 

Using  Fermat's  Principle,  the  reflection  point  will  be  determined 
by  that  v  which  minimizes  s(v) 

s(v)  =  | I ( v) j  +  | d(v)|  (53) 

where 

,  2  1/2 

I  I(v)|  =  [(acosv-xsr  +  (bsinv-y$;  ]  (54) 

1  /2 

|d(v)|  =  [(xc-acosv)2  +  (yc-bsinv)2]  .  l5j 

The  extrema  of  s(v)  can  be  found  by  computing  the  roots  of 
s'(v).  Since  s(v)  is  a  non-linear  function  it  will  be  necessary 
to  use  a  numerical  procedure,  and  Newton's  method  will  work  well 
here.  Substituting 


a  =  cosv 
B  =  sinv 


(56a, b) 


(57) 


h(v)  =  (aorxc)2  +  (bg-yc)2 


g(v)  =  (ac rx$)2  +  (bg-ys)2  (58) 

then 

s(v)  =J~h  +  Jg~  (59) 

and 

s'(»)  =  iiiU-tiLal  (60, 

2  hg 

s"(v)  J2h;'Mh!l2l93/2  *  j29"9-(9’)W/Z  .  ,6)) 

4(hg)J/2 

Now  using  Newton's  method 
s'(v) 

\j  ~  \J  —  - -J - 7- 

n+l  n  s"(v) 

by  observation  of  the  physical  problem  it  is  apparent  that  s'(v) 
will  have  at  most  four  and  at  least  two  roots.  The  four  root  case 
(two  maximums  and  two  minimums)  occurs  when  the  source  and  observa¬ 
tion  points  are  near  enough  so  that  the  path  through  the  cylinder 
is  shorter  than  the  path  to  the  ends  (Figure  15).  In  any  case  the 
second  minimum  will  always  lie  outside  the  visible  region. 

Newton's  method  requires  that  an  initial  value  of  v  be  found 
to  start  the  iteration  process.  If  the  initial  value  is  close  to 
the  desired  result  then  successive  iterations  will  always  converge 
to  the  desired  root  and  the  other  roots  need  not  be  calculated. 

No  method  has  been  determined  to  pick  such  an  initial  value.  The 
other  alternative  is  to  choose  a  result  which  is  known  to  converge 
to  the  desired  root  most  of  the  time.  The  value  then  obtained  is 
checked  to  insure  that  the  true  reflection  point  has  been  determined 
if  not,  then  a  slower  but  more  reliable  method  can  be  employed,  or 
a  new  initial  value  can  be  chosen  based  on  some  different  criterion. 

Numerical  analysis  has  shown  that  the  value  of  v,  given  in 
Equation  (63)  will  converge  to  the  proper  root  in  almost  all  cases. 

A  wide  variety  of  cylinder  source  and  observation  points  were  tested 
and  this  criterion  was  found  to  fail  only  for  a  certain  orientation 

when  ~  »1  or  when  one  point  was  in  the  far  field 


Figure  15.  Various  examples  of  the  extremia  for  the 
cylinder,  observer  ray  path. 

source, 

-  ‘-1  R) 

(62a) 

•«  ■  R) 

(62b) 

v +v„ 
s  c 

Kvcl  >  11 

(63a) 

vs+vc 

2  +ir 

Kvcl  > 

(63b) 

Once  a  root  has  been  found  using  v,  it  is  necessary  to  deter¬ 
mine  if  it  is  the  true  reflection  point.  There  are  two  possibilities 
for  error:  1)  it  has  converged  to  a  maximum,  2)  it  has  converged 
to  a  minimum  not  In  the  visible  region.  For  the  first  case  it  is 
necessary  only  to  check  s"(v)>0.  In  the  second  case,  it  is  necessary 
to  determine  whether  the  lines  (x  y  )  to  (x  y  )  or  (x  y£)  to  (xryr) 
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Since 


intersect  the  cylinder  (except  at  the  point  (xryr). 


x  =  A  cos  v 
r 

(64) 

and 

yr  =  B  sin  v 

(65) 

then  from  formula  (49) 

(m2  +  x2  +  2mcx  +  (c2-b2)  =  0 

3 


(66) 


Since  one  root  of  the  equation  is  already  known  to  be  x,=x  ,  then 
by  division  the  second  root  is 

?  h2 
2mc  -  x  (m  +  -2*.) 
r  2' 

x2  = - ,  ,2  ■  '  (67> 


m2  + 


b£ 

T 

a 


The  intersection  points  have  now  been  calculated  and  it  is  a  simple 

matter  to  test  whether  or  not  the  ray  path  intersects  the  cylinder 

testing  whether  x  <  x,  <  x  or  x  >  x„  >  x  . 

r  2  s  r  2  s 

D.  Tangent  Points 


For  the  case  where  the  source  is  hidden  by  the  cylinder,  the 
reflection  point  does  not  exist.  In  that  case  it  is  necessary  to 
compute  the  tangent  points  on  the  cylinder  since  that  is  where  the 
creeping  waves  will  couple  on  and  shed  from  the  surface.  These 
points  (x,y)  are  given  by  [2] 


12 


Je2-2-2- 2  -2.2 


a2b2x$  ?  z2ysJa‘-y^bts^-atbt 

jrj 


(a‘v;  +  b2x2) 


(68) 


aVys  i  b2x$. 


J.2..2 


a“y^+b2x2-a2b2 


( 6h) 
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the  tangent  angles  can  then  be  found  from 


1,2 


tan 


(70) 


once  the  tangent  points  are  known  in  the  (x,y)  plane,  the  clockwise 
and  counter  clockwise  paths  around  the  cylinder  can  guickly  he  de¬ 
termined.  The  only  unknowns  at  this  point  are  the  /.  coordinates 
of  the  coupling  points.  The  z  coordinate  can  be  determined  by  making 
use  of  the  fact  that  the  elliptic  cylinder  is  a  developed  surface. 


+  ,  (zc"zs) 

zs  "  dos  ( d  +t+d  ) 
„  „  os  oc' 
c  oc 


(71) 


where 


"os  ■J(Vxi/*  <Vyi/ 

" oc  SJ(X <yc-yif 

Q, 


t  =  /  v/a?sin2(v)  +  b2cos2(v) 


dv 


Results 


(72a) 

(72b) 

(73) 


Figures  17,  18  show  plots  of  the  coupling  between  two  X/2  di¬ 
poles  using  the  geometry  given  in  Figure  16.  The  data  is  plotted 
with  the  field  computed  by  the  NEC  Basic  Scattering  Code  [5,12], 

The  NEC  is  a  computer  program  developed  at  the  Ohio  State  Uni 
versity  ElectroScience  Laboratory.  It  computes  the  far  field  scat¬ 
tering  of  electric  and  magnetic  sources  in  the  vicinity  of  elliptic 
cylinders  and/or  finite  plates.  This  code  was  chosen  as  a  basis  for 
comparison  because  of  its  flexibility  and  proven  accuracy.  Since 
the  fields  computed  by  NEC  are  scattered  electric  fields,  they  can 
he  used  for  comparison  only  with  the  following  points  in  mind. 
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Figure  16.  3-D  scattering  configuration  involving  an  infinitely 

long  elliptic  cylinder. 


1)  The  test  dipole  must  be  oriented  to  pick  up  only  the 

components  of  the  field  computed  by  NEC  (F  or  E  ). 

T  t 

2)  the  test  dipole  must  be  far  enough  from  the  cylinder 
and  the  source  dipole  that  the  far  field  approximations 
apply. 

3)  Since  the  hybrid  code  is  computing  coupling,  only  the 
normalized  pattern  can  be  compared.  This  is  possible 
because  the  test  dipole  is  in  the  far  field  and  short 
compared  with  its  distance  from  the  cylinder.  Thus 
the  incident  field  can  be  considered  constant  across 
its  length.  With  a  constant  aggument  in  the  integrand 
the  incident  field  and  coupling  differ  only  by  a  con¬ 
stant  value. 

It  was  also  desired  to  show  that  the  hybrid  method  works  well 
when  both  dipoles  were  near  the  cylinder. 


A  Green's  function  solution  for  the  coupling  between  two  axial 
PWS  dipoles  in  the  presence  of  an  elliptic  cylinder  was  provided 
by  Ersoy  [13].  This  procedure  is  very  similar  to  the  hybrid  technique 
described  in  this  paper.  A  delta  impedance  matrix  representing  the 
cylinder  effect  is  found  via  a  moment  method  procedure.  The  method 
incorporates  the  cylindrical  Green's  function  in  the  kernel  of  the 
integral  equation.  These  eigenfunctions  account  for  the  presence 
of  the  cylinder.  Since  the  eigenfunctions  are  in  the  form  of  a  sum¬ 
mation,  accurate  results  are  obtained  when  sufficient  terms  are  in¬ 
cluded.  This  will  be  seen  to  pose  some  problems. 


Figures  19  and  20  show  how  the  hybrid  and  the  eigenfunction 
solution  compare  using  the  geometry  of  Figure  16.  Agreement  is  very 
good  except  in  the  deep  shadow  region.  Although  the  eigenvalue  solu¬ 
tion  is  analytically  exact,  it  could  suffer  from  the  aforementioned 
numerical  problems.  Thus  another  method  of  checking  the  result  was 
desired. 


F.  Experimental  Data 

The  only  other  method  immediately  available  for  comparison 
was  the  measurement  of  actual  antenna  systems.  To  avoid  the  problems 
associated  with  balanced  and  unbalanced  antenna  feeds,  the  model 
was  constructed  on  a  ground  plane  and  monopoles  were  used  instead 
of  dipoles.  Figure  21  shows  the  model  used.  A  five  inch  diameter 
cylinder  was  mounted  on  the  ground  plane,  and  an  aluminum  "skirt" 
was  put  on  the  plane  to  reduce  edge  diffractions.  The  monopoles 
were  constructed  from  1/8"  diameter  brass  rod  and  cut  to  2.5"  in 
length.  The  monopoles  were  A/4  long  at  1.1811  GHz  Ekleman  p]used 
a  similar  arrangement  to  measure  input  impedance.  He  showed  that 
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HYBRID  SOLUTION 
EIGENFUNCTION  SOLUTION 


since  the  monopoles  are  moderately  thick,  it  was  necessary  to  lengthen 
the  monopoles  in  the  computational  model  by  1/16"  to  account  for 
tip  effects.  This  extra  length  was  also  included  here. 


The  test  equipment  used  was  designed  to  accurately  measure 
the  transmission  coefficient  of  a  two  port  network.  The  two  antenna 
system  can  be  viewed  as  such  a  network  with  the  input  port  being 
the  terminals  of  the  transmitting  antenna  and  the  output  port  being 
the  terminals  of  the  receiving  antenna.  The  coupling  Z,?  is  related 
to  the  transmission  coefficient  T^  by  1 


T 


12 


(74) 


where  Z ^  is  the  input  impedance  of  the  receiving  antenna. 


Data  was  taken  over  frequencies  from  1.  to  1.8  GHz.  In  terms 
of  wavelengths  this  provided  a  variety  of  antenna  lengths  and  cyl¬ 
inder  sizes.  The  measured  and  computed  values  are  shown  in  Figures 
23  and  24.  Agreement  between  the  two  sets  of  data  is  good,  with 
the  measured  points  showing  good  spread  around  the  computed  curves. 
Since  the  losses  were  so  great(|T,2|  ~  44  dB),  it  was  necessary  to 
use  very  high  gain  settings  of  the^test  equipment  amplifier.  This 
caused  the  system  to  be  very  sensitive  to  outside  influences  such 
as  other  objects  in  the  anechoic  chamber.  In  spite  of  this  sensi¬ 
tivity,  several  sets  of  data  were  taken  and  good  repeatability  of 
the  data  was  observed.  Thus  the  hybrid  method  has  been  shown  to 
accurately  predict  the  coupling  between  antennas  near  a  cylinder. 


Figures  25,  26,  and  27  show  the  coupling  between  various  ori¬ 
entations  of  dipoles  near  infinite  cylinders.  Note  that  in  Figure 
25  the  free  space  and  diffracted  curves  differ  only  by  a  constant 
in  the  far  field.  This  occurs  because  the  path  around  the  cylinder 
does  not  vary  in  that  region.  Note  also  that  in  Figures  25  and  26 
a  peak  is  seen  in  the  coupling.  This  occurs  because  in  the  extreme 
near  field  the  field  is  highly  attenuated  due  to  the  long  path  around 
the  cylinder.  In  Figure  27  this  peak  does  not  occur  because  the 
path  around  the  cylinder  is  always  getting  longer. 
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Figure  22.  Block  diagram  of  experimental  measurement  setup. 
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Figure  23.  Magnitude  of  the  transmission  coefficient  for  the 
monopoles  used  in  the  experimental  setup. 
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CHAPTER  VI 

MONOPOLES  ON  CURVED  SURFACES 


The  preceeding  section  presented  a  general  method  for  comput¬ 
ing  the  currents  on  thin  wire  antennas  near  a  cylinder.  A  related 
problem  of  great  interest  is  when  the  antenna  is  attached  to  the 
cylinder  and  driven  against  it.  The  simplest  example  of  this  is 
the  case  of  a  monopole  which  is  normal  to  the  surface  at  the  attach¬ 
ment  point. 

It  would  be  desirable  to  continue  to  employ  the  piecewise  sinu¬ 
soidal  Galerkin  method  since  the  solution  would  be  compatible  with 
the  hybrid  method  already  presented.  To  accomplish  this  some  dif¬ 
ficulties  must  first  be  resolved. 


The  piecewise-sinusoidal  functions  are  efficient  because  they 
are  very  similar  to  the  actual  current  distribution  on  a  wire.  One 
of  the  most  noteable  features  is  that  the  current  goes  to  zero  (as 
expected)  at  the  ends  of  the  wire.  When  this  wire  is  attached  to 
a  surface  this  is  no  longer  true.  Since  the  current  is  in  general 
not  zero  at  the  feed  point  of  a  monopole  special  techniques  will 
have  to  be  developed  in  order  to  handle  this  problem. 


A.  Image  Theory  Solution 


If  the  radius  of  curvature  of  the  surface  is  large,  then  image 
theory  can  be  used  to  compute  the  currents  without  difficulty.  This 
case  is  shown  in  Figure  27.  In  this  figure  the  wires  have  been 
broken  up  into  segments  and  the  corresponding  PWS  modes  have  been 
included  and  numbered  for  reference.  Writing  the  equations  of  Figure 
27  gives 


Z11I1  +  Z12 1 2  +  Z13 1 3  +  Z14!4  +  Z15I5  +  Z16I6  +  Z17J7 
Z21 1 1  +  Z22 1 2  •' 


Z7 1 1 1  +  Z72*2 


Z  77 1 7 


(75) 
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Figure  27.  Example  of  the  segments  and  mode  structure  on  a 
monopole  and  its  dipole  image  equivalent. 


from  symmetry  it  is  known  that 
I 


‘1 


^7»  ^  *6’  *3  =  *5 


and  thus  in  matrix  form 


[V]  = 


=  < 


Z11 

Z12 

Z13 

Z14_ 

_Z17 

Z16 

Z15 

O' 

Z21 

Z22 

Z23 

Z24 

Z27 

Z26 

Z25 

0 

31 

Z32 

Z33 

Z34 

+ 

Z37 

Z36 

Z35 

0 

Z41 

Z42 

Z43 

Z44 

_Z47 

Z46 

Z45 

0 

[I] 


(76) 


(77) 
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or 


! Vi  -  <  rzpi  +  [az]}  m 


'781 


Several  important  observations  can  be  drawn  from  (781.  First, 
the  Zp  matrix  is  merely  the  standard  impedance  matrix  which  would 
be  constructed  if  the  antenna  had  been  defined  using  only  segments 
1-5  in  free  space.  Secondly,  the  [AZ]  matrix  consists  of  couplings 
only  between  modes  on  opposite  sides  of  the  ground  plane.  These 
couplings  are  the  same  as  the  coupling  between  modes  in  [Zp]  if  only 
reflected  ray  paths  were  considered  since  the  modes  below  the  ground 
plane  are  merely  images  of  the  segments  above. 

Thus  instead  of  formulating  the  problem  using  image  theory, 
the  problem  can  be  solved  directly  by  computing  a  free  space  imped¬ 
ance  matrix  and  a  perturbation  impedance  matrix  using  reflection 
theory.  This  result  is  not  too  surprising  when  one  considers  the 
following.  Each  term  Z  of  the  matrix  [Z]  is  obtained  by  computing 
the  coupling  between  a  test  dipole  and  an  expansion  dipole  with  no 
regard  to  the  rest  of  the  antenna.  Thus  except  for  the  case  where 
the  test  or  expansion  dipoles  is  in  contact  with  the  surface,  the 
theory  presented  in  Chapter  IV  would  suffice.  The  discussion  on 
image  theory  was  necessary  only  to  justify  the  existence  of  the  mode 
which  "straddles"  the  surface  (mode  4)  and  to  provide  a  means  of 
including  it  in  the  current  theory.  The  ground  plane  can  be 
replaced  by  a  curved  surface  as  follows.  Recall  that  each  term 
(Z  )  of  the  [Z1  matrix  represents  the  total  coupling  between  dipole 
modes  m  and  n,  regardless  of  the  structure  of  the  rest  of  the  an¬ 
tenna.  Thus  for  modes  near  a  curved  surface  this  includes  rays 
diffracted  or  reflected  by  the  surface.  If  the  radius  of  curvature 
of  the  surface  is  not  large,  then  it  is  necessary  to  use  GTD  or  GO 
in  order  to  obtain  an  accurate  result.  The  one  exception  to  this 
is  when  one  of  the  wire  segments  is  quite  near  the  surface  and  ex¬ 
tremely  short,  then  the  surface  still  appears  as  a  ground  plane  and 
image  theory  can  be  used.  Thus  in  the  case  of  a  monopole,  the  [ zJ 
matrix  can  be  constructed  using  GO  or  GTD  for  all  terms  except  for 
the  mode  which  "straddles"  the  surface.  This  "straddling"  mode  can 
be  adequately  handled  with  image  theory,  with  the  restriction  that 
it  must  be  much  shorter  than  the  radius  of  the  cylinder. 

B.  Numerical  Implementation 

Several  problems  were  encountered  when  the  above  technique 
was  used  in  a  computer  program.  Consider  the  coupling  between  modes 
a  and  b  which  lie  on  a  wire  normal  to  the  surface.  In  this  case 
the  reflection  point  is  where  a  straight  line  extending  the  dipole 
to  the  cylinder  would  touch.  Since  GO  or  GTD  only  includes  fields 
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transverse  to  the  direction  of  propagation,  the  Z  .  term  would  be 
zero.  This  occurs  because  the  field  must  be  dotted  with  the  direc¬ 
tion  of  the  segment.  Since  the  direction  of  the  segment  is  along 
the  ray  path,  this  dot  product  is  zero.  In  practice,  this  problem 
is  resolved  by  putting  the  test  dipole  on  the  surface  of  the  wire 
and  the  expansion  dipole  along  the  center  of  the  wire.  This  however 
raises  two  other  questions. 

First,  since  the  wire  is  assumed  to  be  round,  the  test  dipole 
can  be  placed  at  any  position  on  the  outer  surface  without  affecting 
the  geometry  with  respect  to  the  wire.  However  this  affects  the 
geometry  of  the  two  dipoles  with  respect  to  the  curved  surface. 

In  fact  the  only  cases  where  this  placement  does  not  have  any  effect 
is  when  the  surface  at  the  reflection  point  is  spherical  or  planar. 
Numerical  experiments  with  a  thin  wire  normal  to  an  elliptic  cylinder 
were  performed  using  the  geometry  in  Figure  28.  These  results  (Fig¬ 
ures  29,  30)  showed  that  for  cylinders  of  reasonable  size  the  slight 
change  in  geometry  does  affect  the  computed  coupling,  this  effect 
was  quite  small.  Thus  the  placement  of  the  test  dipole  was  arbi¬ 
trarily  chosen  so  that  0  =  0°. 

Another  possible  source  of  numerical  error  arises  from  the 
radial  wire.  One  of  the  basic  assumptions  used  to  compute  the  cou¬ 
pling  in  this  paper  is  that  the  wires  are  electrically  thin.  This 
is  not  a  limitation  of  the  moment  method  but  a  requirement  imposed 
to  simplify  the  computational  problems.  Since  the  wires  are  on  the 
order  of  .006A  or  less,  the  angle  of  incidence  is  quite 
small  even  for  wires  which  are  close  to  the  surface.  This  would 
appear  to  pose  great  numerical  problems.  Again  numerical  experiments 
show  that  this  generally  does  not  lead  to  incorrect  results.  The 
reason  for  this  is  as  follows.  It  is  well  known  that  the  endfire 
field  of  a  dipole  is  only  significant  in  the  near-field.  Examination 
shows  that  as  G  +  0,  E  -*■  0  and  E  falls  off  as  1/R  ,  thus  the  field 
quickly  dies  away.  Numerical  resalts  showed  that  for  wire  radii 
as  small  as  .0005A  with  6  digit  precision  the  field  became  negligibly 
small  before  numerical  problems  due  to  0  *  0  were  encountered. 

Another  problem  resulted  from  the  tip  of  the  mode  which  touched 
the  surface,  mode  3  in  Figure  27.  When  the  self-coupling  of  mode 
3  is  being  computed,  the  distance  from  the  source  point  (on  the  out¬ 
side  of  the  wire)  and  the  observation  point  (on  the  axis  of  the  wire) 
is  quite  small.  Since  the  source  point  is  a  caustic,  the  field 
varies  rapidly  in  this  region.  This  rapidly  changing  field  requires 
that  the  expansion  mode  (along  the  axis)  be  broken  up  into  many 
integration  points  for  the  numerical  integration  to  yield  an  accurate 
result. 
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TOP  VIEW 


I  igure  28.  Location  of  the  test  and  expansion  dipoles  on  a  thin 
wire  mounted  on  a  cylinder. 
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Figure  29.  Real  input  impedance  of  a  monopole  as  a  function  of 
the  location  (9)  of  the  test  dipole. 


c. 


Results 


Figures  31,32  show  the  computed  input  impedance  of  a  monopole 
on  a  circular  cylinder.  The  length  was  chosen  such  that  the  monopole 
is  resonant  when  the  cylinder  becomes  a  flat  ground  plane.  Figure 
33  shows  how  the  hybrid  technique  compares  with  experimental  data 
supplied  by  Wang  [14].  Note  that  except  for  the  classical  resonance 
shift  the  curves  "track"  well. 

While  the  image  theory  solution  of  the  monopole  problem  gives 
accurate  results,  it  has  the  drawback  of  requiring  that  one  mode 
"straddle"  the  surface.  This  mode  also  has  a  length  restriction 
which  usually  increases  the  number  of  required  segments.  Finally 
it  is  a  non-physical  situation  and  thus  is  difficult  to  work  with 
on  an  intuitive  basis.  These  difficulties  could  hinder  future  work. 
Because  of  these  problems  an  alternate  technique  referred  to  as  the 
half  mode  solution  was  developed  which  alleviated  the  above  problems. 

D.  Half  Mode  Solution 

This  solution  sets  up  only  half  of  the  mode  that  originally 
"straddled"  the  surface.  The  couplings  computed  as  a  result  of  this 
half  mode  are  not  valid  until  the  [Zrl  and  [AZ1  impedance  matrices 
are  added  together  to  produce  the  total  impedance  matrix  [Z]. 

As  stated  earlier,  the  coupling  between  two  PWS  dipoles  is 
actually  computed  via  the  coupling  between  four  PWS  monopoles.  These 
couplings  are  indicated  as  arrows  in  Figure  34.  The  point  charge 
(due  to  the  non-zero  current  at  one  end)  that  would  result  from  each 
monopole  is  neglected  since  it  will  be  cancelled  out  when  these 
monopoles  are  connected  as  dipoles.  Now  consider  the  coupling  be¬ 
tween  a  PWS  dipole  and  a  PWS  monopole  on  a  ground  plane  in  Figure 
35  where  the  image  theory  couplings  have  been  shown  as  dotted  lines. 
Observe  that  via  image  theory,  this  is  identical  to  the  coupling 
between  two  dipoles,  but  without  explicitly  defining  the  image  arm 
of  the  monopole.  This  forms  the  basis  by  which  the  half  mode  solu¬ 
tion  works.  Careful  examination  shows  that  this  is  just  a  logical 
extension  of  the  image  theory  solution.  Thus  it  was  not  necessary 
to  define  all  of  the  mode  that  straddled  the  surface  but  only  the 
portion  above.  The  only  modification  being,  any  coupling  between 
half-modes  (i.e.,  monopoles  on  the  surface)  must  be  doubled  to  con¬ 
vert  them  to  dipole-dipole  couplings. 

The  two  impedance  matrices  generated  by  this  method  have  a 
peculiar  quality  which  can  be  explained  as  follows.  Letting  the 
half  mode  of  Figure  35  be  labeled  x  and  the  dipole  mode  be  labeled 
y,  the  following  observation  can  be  made.  The  terms  z  and  z 
of  the  matrix  [Zpl  are  not  valid  by  themselves  (i.e.,  Without  adding 
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Figure  34.  Graphical  representation  of  the  monopole-monopole 
couplings  required  for  a  dipole-dipole  coupling. 


Figure  35.  Image  equivalent  of  Figure  34. 
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Az  or  Az  ).  This  arises  because  each  term  in  [Zr]  is  by  definition 
a  flTpole  tcrdipole  impedance  but  z  and  z  are  monopole-dipole 
impedances  (without  the  point  charge  contribution) .  However  since 
Azxy  and  az  contain  the  corresponding  monopole-dipole  impedances 
(again  without  point  charge  contribution)  their  sum  does  define  a 
dipole-dipole  impedance  this  giving  v  valid  impedance  matrix.  The 
computer  programs  used  to  generate  the  half  mode  solution  were  changed 
to  construct  a  solution  using  the  half-mode  technique.  The  results 
were  identical  (within  numerical  error)  to  those  produced  by  the 
image-theory  solution.  The  drawback  to  this  method  is  that  it  is 
not  directly  compatible  with  Richmond's  thin  wire  programs.  These 
routines  were  never  intended  to  be  used  on  a  system  where  such  a 
mode  could  exist  and  thus  they  would  have  to  be  changed.  For  the 
purposes  of  this  report,  the  half  mode  solution  was  only  implemented 
well  enough  to  run  a  few  cases  and  demonstrate  that  it  would  work. 

The  solution  was  presented  only  as  a  stepping  stone  to  future  work 
in  this  field. 


CHAPTER  VII 

SUMMARY  AND  CONCLUSION 


The  hybrid  technique  of  Ekelman  has  been  expanded  to  handle 
a  more  general  set  of  thin  wire  antennas  and  scatterers  than  was 
previously  possible.  The  uniform  GTD  of  Pathak  et  al.  [3  Iwas  used 
to  provide  a  means  of  analyzing  systems  where  antennas  are  arbitrarily 
oriented  near  an  infinite,  perfectly-conducting  elliptic  cylinder, 
and  a  new  hybrid  image  theory  technique  was  developed  to  analyze 
monopoles  mounted  on  the  cylinder.  Both  of  these  new  solutions  re¬ 
tain  the  format  and  conventions  set  up  by  Richmond  and  are  completely 
compatible  with  his  powerful  moment-method  computer  code. 

The  free  space  impedance  matrix  was  computed  using  the  con¬ 
ventional  moment  method  technique.  A  perturbation  matrix,  called 
the  delta  impedance  matrix  after  Thiele  et  al.  was  then  added  to 
the  free  space  matrix.  This  delta  matrix,  contained  the  coupling 
between  elements  due  to  GTD  type  rays  which  were  reflected  or  dif¬ 
fracted  by  the  cylinder.  Thus  the  effects  of  the  cylinder  were 
accounted  for  in  the  final  solution. 

The  first  problem  solved  was  that  of  the  piecewise  sinusoidal 
dipoles  in  the  presence  of  the  aforementioned  cylinder.  A  dipole 
with  a  single  PWS  current  distribution  is  the  single  most  important 
type  of  antenna  that  can  be  considered  since  each  element  of  the 
impedance  matrix  is  merely  the  coupling  between  two  such  dipoles. 

The  solution  obtained  was  shown  to  be  accurate  by  comparing  it  with 
three  independent  techniques.  The  NEC  Basic  Scattering  Code  developed 
at  The  Ohio  State  University  was  used  to  check  general  orientations 
in  the  far  field  results.  A  moment  method  eigenfunction  technique 
was  used  to  check  the  results  when  both  dipoles  were  in  the  near 
field.  The  eigenfunction  solution  had  problems  in  the  deep  shadow 
region  so  measured  data  was  used  to  check  that  case.  All  of  the 
independent  results  showed  the  hybrid  method  to  be  accurate  and  re¬ 
liable. 

The  second  class  of  problem  solved  was  a  monopole  mounted  on 
a  curved  surface.  The  solution  was  obtained  by  image  theory  with 
perturbations  of  the  images  by  the  curved  surface.  This  solution 
is  very  general  and  is  applicable  to  any  curved  or  flat  surface. 

To  demonstrate  its  usefulness,  the  monopole  was  put  on  an  infinite. 
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perfectly  conducting  elliptic  cylinder  and  its  input  impedance  was 
obtained.  This  computed  impedance  was  shown  to  agree  favorably  with 
measured  data.  A  second  method  called  the  half-mode  solution  was 
also  presented  for  solving  the  monopole  problem.  This  method  used 
the  same  theory  as  the  original  but  in  a  manner  which  could  make 
future  work  easier.  While  the  half-mode  solution  was  explained 
fully,  it  was  intended  only  to  present  a  possible  alternative. 

The  implementation  was  carried  only  far  enough  to  demonstrate  that 
it  will  work.  More  work  will  need  to  be  done  if  the  full  potential 
of  this  method  is  to  be  realized. 
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